Structure and static response of small silver clusters to an external electric field 
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The static response properties and the structural stability of silver clusters in the size range 
1 < n < 23 have been studied using a linear combination of atomic Gaussian-type orbitals within the 
density functional theory in the finite field approach. The Kohn-Sham equations have been solved 
in conjuction with a generalized gradient approximation (GGA) exchange-correlation functional. 
A proof that the finite basis set GGA calculation holds the Hellmann-Feynman theorem is also 
included in the Appendix. The calculated polarizabilities of silver clusters are compared with the 
experimental measurements and the jellium model in the spillout approximation. Despite the fact 
that the calculated polarizabilities are in good agreement with both of them, we have found that the 
polarizability appears to be strongly correlated to the cluster shape and the highest occupied-lowest 
unoccupied molecular-orbital gap. 

PACS numbers: 36.40.Cg, 32.10.Dk, 31.15.Ew, 61.46.Bc 



I. INTRODUCTION 

Transition-metal clusters play a dominant role in clus- 
ter physics [l[ . They have attracted the interest of many 
researchers and consequently the number of publications 
on that topic have experienced a dramatic increase over 
the past thirty years because metal clusters exhibit in- 
creasingly interesting structural, electronic, catalytic, as 
well as optical properties [2|, |3( . Likewise, the biomedi- 
cal applications of nanoclusters have experienced a great 
impact in the biomedicine community [J]. 

Among the aforementioned applications and properties 
of metal clusters, their optical properties that result from 
the light-matter interaction at the nanoscale level are an 
area of great current interest [5| . Many of the advances in 
this area could not have been made possible without the 
development of the optical spectroscopy techniques that 
have been indispensable for elucidating the electronic 
structure of clusters. These experimental techniques can 
be divided into two groups, that is, nondestructive and 
destructive methods. In the former methods, also called 
linear response methods, a weak electromagnetic field in- 
teract with the cluster and it absorbs or scatters light 
without undergoing ionization or dissociation whereas in 
the destructive methods the ionization is achieved. The 
nondestructive methods in connection with the linear- 
response theory have been extensively used to calculate 
photoabsorption cross sections and specially static dipolc 
polarizabilities. 

The static dipole polarizability is a physical observ- 
able of metal clusters that has been shown to be closely 
related to the shape and structural geometry. For ex- 
ample, electronic structure calculations of small Si clus- 
ters show that the polarizability is strongly correlated 
with the shape of the clusters @. Likewise, the interplay 
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between theory and experiment is a powerful tool that 
serves to identify which cluster is observed in the ex- 
periments throughout the comparison of the calculated 
polarizabilities with the experimental ones • Moreover, 
the static dipole polarizability is also well-known that it 
is intimately related to the shell electronic structure. For 
example, the noble-metal clusters whose optical proper- 
ties have been extensively studied in literature (mainly 
from the experimental side) present lower static polar- 
izabilities than alkali metals because they are excellent 
examples of spherical shell structure. However, in the 
case of the silver clusters, the influence of d electrons 
in the static polarizability has been less studied at least 
from the theoretical point of view and it deserves more 
investigation. This point, among others, is addressed in 
this work. 

In this article we have employed a density functional 
theory-based (DFT) calculation within the generalized 
gradient approximation (GGA) to properly account the 
strong correlation effect of the localized d electrons and 
charge density inhomogeneities. We have studied the 
structural stability and the static response properties of 
small silver clusters ranging in size from n=l up to n=23, 
where n is the number of atoms forming the cluster. The 
calculated static polarizabilities are only compared with 
the available experimental data since that unfortunately, 
no ab initio quantum-molecular calculations have been 
done so far for the static polarizabilities of silver clus- 
ters in the size range covered by our investigation. The 
available theoretical data are reported in Ref. Q for very 
small silver clusters with 1 < n < 8. The agreement 
between our values and those of Idrobo et al, who com- 
puted the static polarizabilities within the framework of 
the real space finite-different ab initio pseudopotential 
method, is excellent in the case of the GGA approxi- 
mation, which makes more valuable the results obtained 
with the present DFT method. In addition, we have 
compared our results with the jellium model in the spill- 
out approximation and the deviation of our results from 
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this model is explained in terms of the electronic struc- 
ture parameters like highest occupied and lowest unoc- 
cupied molecular orbital (HOMO-LUMO) gaps or struc- 
tural symmetry. We also show and discuss how the d 
electrons affect the static polarizabilities since that when 
n is sufficiently large, i.e. greater than 18 the polariz- 
ability tends to be constant in contrast with the one of 
the alkali metals. The rest of the paper is organized as 
follows. In Sec.|TTl we present the theoretical background 
and the computational details along with the compu- 
tational parameters used in this article. Moreover, the 
structural stability of the silver clusters is studied in de- 
tail in Sec. IIIIl The results of our calculations and the 
influence of the electronic structure observables in the 
static polarizabilities are given in Sec. IIVI We conclude 
with a brief summary of the reported results provided 
by our ab initio calculations in Sec. |Vj In the Appendix 
we prove that the Hellmann-Feynman theorem is satis- 
fied by the finite basis set density functional framework 
in the generalized gradient approximation. 



II. METHOD AND COMPUTATIONAL 
DETAILS 

Traditionally there has been two different ways of com- 
puting the polarizability of molecular systems. Thus, 
the polarizability is identified either as the second-order 
term in the perturbation expansion of the electronic en- 
ergy with respect to the applied external uniform elec- 
tric field or as the linear response of the dipole mo- 
ment to that electric field. Both definitions are equiv- 
alent when the Hellmann-Feynman theorem is satisfied. 
The Hellmann-Feynman theorem holds for an exact so- 
lution of the Schrodinger equation and also for some ap- 
proximate solutions. In particular, it is satisfied by the 
fully self-consistent finite basis set solutions when the 
exchange-correlation (XC) energy is approximated by the 
GGA as shown in the Appendix. 

We have adopted the dipole moment-based definition 
as our working definition because in this expression the 
field occurs only in the first power for the calculation 
of the polarizability instead of the second as in the en- 
ergy expansion. The polarizabilities ({i,j}=x,y,z or 
alternatively |Lj}=l,2,3) are calculated by the finite field 
(FF) method [!| which consists of computing the electric 
dipole moment ^ of a system under the influence of an 
external electric field Fi according to the following finite- 
difference relation 



= km tiM^itljl. (i) 

F=0 2F i 



In the FF method, one of the most crucial problems to 
evaluate the derivatives is the choice of an appropriate 
field strength. Several works have assessed the numeri- 
cal accuracy of polarizability against different field values 
and concluded that the best region of linear response is 



[id , lllj | . For that reason we have used a field strength of 
5 x 10~ 4 a.u. that is applied along the molecular axis. 
At least 7 self-consistent field (SCF) runs, with the field 
strengths and ±Fj, are necessary to calculate the polar- 
izability. Once the polarizability tensor components are 
computed, the mean static polarizability is calculated as 
a = (E f =1 a^)/3 and the polarizability anisotropy is de- 
fined as 
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eneral axis frame or without the second addended 
(6 EfjLi 2-i<j a ij) m t ne coordinate system which makes 
the second-rank polarizability tensor a diagonal. 

With the aim of studying the static response proper- 
ties of small silver clusters, Ag n (2<n<23), we have per- 
formed density functional theory-based calculations con- 
sisting of a linear combination of Gaussian-type-orbitals 
(LCGTO) Kohn-Sham density-functional methodology 
as it is implemented in DEMON-ks3p5 program [12] . All- 
electron spin-unrestricted calculations were carried out at 
the GGA level to take the XC effects into account [l3| . 
Local-density approximation (LDA) sometimes yields in- 
accurate bond lengths and total energies due to the insuf- 
ficiency in describing the strong correlation effects of the 
localized d electrons and charge density inhomogeneities. 
In these regards, the GGA should be a choice better than 



LDA [14J. For this reason, at the beginning of this work 



and to satisfy ourselves that the numerical procedure is 
reliable, we initiate a search of the functional that better 
fitted the calculated bond length of the silver dimer to the 
experimental one. The functional developed by Perdew 
and Wang [l3| gave us a bond length of 2.534 A, that 
is in excellent agreement with the experimental measure 
(2.53350 A) reported in Ref. [Tj|. An orbital basis set 
of contraction pattern (633321/53211*/531+) was used 
in conjunction with the corresponding (5,5;5,5) auxiliary 
basis set for describing the s, p and d orbitals [l6| . In 
demon-ks3p5, the electron density is expanded in aux- 
iliary basis functions which are introduced to avoid the 
calculation of the N 4 scaling Coulomb repulsion energy, 
where N is the number of the basis functions. The grid 
for numerical evaluation of the XC terms had 128 radial 
shells of points and each shell had 26 angular points. Spu- 
rious one-center contributions to the XC forces, typically 
found in systems with metal-metal bonds when using a 
nonlocal functional, are eliminated in a similar way as 
has been done in Ref. [17| . A wide set of spin multiplici- 
ties ranging from 1 to 11 was checked to ensure that the 
lowest-energy electronic configuration is reached. The ge- 
ometries were fully optimized without symmetry and ge- 
ometry constraints using the Broyden-Fletcher-Goldfarb- 
Shanno algorithm [Tsj | . During the optimization, the con- 
vergence criterion for the norm of the energy gradient 
was fixed to 10 -4 a.u. while it was 10~ 7 a.u. for the 



for field strengths ranging from 10 up to 10 a.u. energy and 10 a.u. for the charge density. The ground 
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state structures and some of the lowest-energy isomers of 
the silver clusters studied in this article are illustrated in 
Fig. [U since that the polarizability is closely related to 
the geometrical shape of the cluster. 



III. STRUCTURES OF SILVER CLUSTERS 

Figure[T]and Table[T]show the structures, relative ener- 
gies, first-neighbor distances and vertical ionization po- 
tentials (IP) of the lowest-energy isomers predicted by 
the simulations for Ag„ clusters with 2 < n < 23. The 
IPs have been calculated as the energy difference between 
the neutral cluster and the cation. We assign labels to 
clusters such as n_m, where n indicates the number of 
atoms and the second number gives the rank in increas- 
ing energy order. A huge sampling of trial geometries 
taken from the literature was evaluated [lji[. While for 
these small clusters, it is nearly impossible to search for 
all possible geometries, the detailed search that we have 
carried out gives us some confidence that the ground state 
structures have been found. 

For clusters with size varying from 2 up to 13, we opti- 
mized geometries and calculated energies for 78 isomers. 
We only reported the lowest-energy structures here, but 
their Cartesian coordinates in Bohr units are available 
upon request to the authors. We focus on isomers lo- 
cated 50 KJ/mol (w 0.5 eV) or less above the lowest- 
energy structure because this value is considered to rep- 
resent an error bar for relative energies computed with 
GGA functionals [2(J. Whenever is possible we compare 
the calculated IPs with the experimental ones to rule 
out the possible candidates to the ground state geome- 
try. We will focus our attention mainly on the lowest 
energy cluster geometries which are different form the 
cluster structures reported in previous experimental or 
theoretical studies. For the rest of clusters, we simply 
just describe the lowest energy structures. In Table [H 
other low energy isomers are also listed for comparison. 

For the clusters with size varying from 14 up to 23, we 
took the initial guess structures from Ref. [2l| ■ Moreover, 
in the early stage of the geometry optimization process 
and with the aim of speeding up the calculations, the 
structure of the silver clusters was first optimized in con- 
junction with a 17-electron scalar relativistic model core 
potential designed for the adequate description of the 
silver dimer bond length [12] ■ Once the geometry of the 
cluster was converged for the model core potential, an all- 
electron structural-relaxation calculation was performed 
leading to the current lowest-energy structures showed in 

Fig.ru 

For the silver trimer, we have studied only one struc- 
ture (isosceles triangle) because the minimum structure 
is well established in the literature [H, US] • In the case 
of Ag4, the lowest energy structure that we have found 
is the planar rhombus. The lowest-energy structures for 
Ag3 and Ag4 predicted by our calculations are in good 
agreement with other ab initio results [l9l . l23l l2~ij . 



The ground state structure found for the pentamer is 
a trapezoid with C2 V symmetry. In literature there is a 
controversy about the assignment of the structural min- 
imum. Some theoretical studi es p redict as the ground 
state the planar structure [l9|, |24|. [25| while other ones 
predict the trigonal bipyramid [23j |, On the experimen- 
tal side, it is worthwhile to mention that different ex- 
perimental ESR and Raman spectra of Ag 5 have been 
interpreted by both trigonal bipyramid [261 ] and planar 
geometries [2Jj , respectively. As shown in Table Q] the 
vertical IP calculated for the trapezoid agrees quite well 
with the experimental result [28| and it is what underpins 
our predicted lowest-energy structure. 

Four structures have been optimized for the Ag 6 clus- 
ter [29l | . The planar trapezoidal D^h structure resulted 
to be the most stable but only 0.079 eV lower in energy 
than the pentagonal pyramid. Both structures exhibit 
IPs very close to the experimental one and therefore, our 
ab-initio calculations do not allow prediction of only one 
geometry as the structural minimum. 

In the case of the silver heptamer, the pentagonal 
bipyramid (D§h symmetry) is predicted as the lowest en- 
ergy structure while the tricapped tetrahedron is 0.17 eV 
higher in energy. Although the relative energy and the 
ionization potentials do not allow a clear prediction about 
the minimum geometry, however most of the first princi- 
ples calculations prior to this article have also obtained 
the penta gon al bipyramid as the fundamental structure 

For the silver octamer, we have decided to optimize as 
a good candidate to the structural minimum the follow- 
ing isomers: a T)2d dodecahedron, which can be viewed as 
a distorted bicapped octahedron, a T^ tetracapped tetra- 
hedron and a C s pentagonal bipyramid. Our calculations 
stabilize first the dodecahedron, secondly but very close 
in energy the T^ structure and finally the pentagonal 
bipyramid with AEdft— : 0.181 eV. The comparison be- 
tween the calculated IPs and the experimental measure- 
ments reported in Ref. [281 ] favor more the T^ structure 
than the T>2d one, however the assignment of the struc- 
tural minimum is clearly inverted if we use the experi- 
mental IP (6.40 eV) measured in Ref. [3l|. Our predicted 
minimum structure is also supported by a very recent 
reference where the authors determine the lowest-energy 
structure of Ags by comparison of the optical spectra 
provided by time-dependent DFT with the experimental 
findings [32|. 

Fifteen structures have been optimized in the case of 
the Agg cluster. A tricapped-distorted octahedron (C s ) 
was found for the lowest-energy structure which is in 
good agreement with the reported structure in Ref. [lj|, 
while in Refs. [23|, [24[ the ground state geometry is the bi- 
capped pentagonal bipyramid. To satisfy ourself that the 
bicapped pentagonal bipyramid is not the structural min- 
imum, we studied a wide range of pentagonal bipyramid 
structures capped with two atoms in different positions. 
After the geometry optimization, the final geometry is 
slightly distorted in most of the cases. The structures 
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are plotted in Fig. Q] and denoted as 9_1, 9_2, 9_3, 9_4, 
9_5, 9_6, 9_8, 9_9, 9_10, 9_11, and 9_12. It is worthwhile 
to mention that from 9_0 up to 9_6, the structures are 
very close in energy and only the IPs of the pentago- 
nal bipyramid structures agrees quite well with the ex- 



perimental ones, so in gas-phase experiments is likely to 
obtain the bicapped pentagonal bipyramid as the lowest- 
energy structure. 

Contrary to what happened in the case of nanomers, 
the lowest-energy geometry of Agio cluster is a pentag- 
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onal bipyramid-shaped structure, namely, a tricapped 
pentagonal bipyramid. With the aim of testing the va- 
lidity of this structure, we have also optimized a very 
distorted pentagonal bipyramid but it finally converged 
to the 1CL0 structure. The second structure is a bi- 
capped square antiprism and it is only 0.080 eV higher in 
energy with respect to the ground state. This structure 
has been predicted as the fundamental one in Ref. [23j . 
nevertheless the comparison of the calculated IPs re- 
ported in Table [J with the experimental value clearly 
show that the 10_0 structure is the best candidate to 
be the minimum structure. 

The guess structures for Agn and Agi2 were obtained 
by adding atoms to the pentagonal bipyramid structure 
or by removing them from a 13-atom cuboctahedron 
or from the icosahedral packing. Thus, the lowest-energy 
structures come from the pentagonal bipyramid shape in 
both cases while the structures coming from the cuboc- 
tahedron or the icosahedral packing are less favored en- 
ergetically. In the case of the Agn geometries, we have 
studied 9 isomers but except for the 11_1 structure, the 
other ones are energetically far from the lowest-energy 
structure by an amount greater than 0.5 eV. While the 
relative energies clearly show that the fundamental struc- 
ture is the one labeled with 11_0 (AE^p ^ 11 - 1 = 0.268 
eV), however the calculated IP (A IP 11 - 1 =0.02 eV) pre- 
dicts the 1 1_1 as the minimum structure, and conse- 
quently these two structures would be probably observed 
in experiments. The same situation occurs for Agi2 but 
in this case more structures can be observed in low- 
temperature experiments. 

Three structures have been selected as possible candi- 
dates for the structural minimum of the thirteen-atom 
silver cluster. They are the icosahedral geometry (1^ 
symmetry), the cuboctahedron cluster (Oh symmetry), 
and the T)2h structure which is a compact portion of the 
bcc crystal lattice capped with four atoms as it is illus- 
trated in Fig. [TJ After the relaxation of the structures, 
the ground state geometry predicted by our calculations 
is the icosahedral structure. Although the calculated IPs 
are in general in a relatively good agreement with the 
experimental ones, however they do not add too much 
information to the determination of the lowest-energy 
structure. The relative energies of the isomers collected 
in Table Q] together with the fact that most of the articles 
devoted to the study of the structural properties of small 
silver clusters predict the icosahedral structure as funda- 
mental one, give us some confidence that the icosahedral 
packing is valid for the Agi 3 cluster [23|, [H, HI] . 

As we comment above, for the rest of the clusters, that 
is, from Agi4 up to Ag23 we took the initial guess struc- 
tures from Ref. [111. They follow an icosahedral growth 
sequence capped with a variable number of atoms de- 
pending on the cluster size. In general, the calculated 
IPs are in good agreement with the experimental mea- 
surements which make more valuable the geometries opti- 
mized with the DEMON-ks3p5 program using as starting 
point the structures predicted in Ref. [2lj . 



IV. RESULTS AND DISCUSSION 

The main results of our theoretical calculations con- 
cerning the static response properties in conjunction with 
some selected electronic structure properties of small sil- 
ver clusters are collected in Table HTl and Fig. [2] Hereafter, 
the reported results are only for the ground state struc- 
tures. We observe that the calculated polarizabilities are 
in good agreement with the experimental measurements 
reported in Ref. [35|. They approach each other as the 
cluster size increases because the experimental data are 
less reliable as the cluster size decreases. The polariz- 
abilities reported by Fedrigo et al. [35j were measured at 
T= 10 K whereas our calculated results are given at T= 
K. Thus, according to the following relation [36j in the 
low electric field limit 



Qfoff = a + 



3K B T 



(3) 



for clusters having a permanent dipole moment /z, the 
effective measured polarizability at temperature T is ex- 
pected to be greater than average polarizability. However 
we have found that the contribution of the second term 
in Eq. (^-calculated at T= 10 K and for the dipole mo- 
ments collected in Table ITU-is negligible and consequently 
the influence of the temperature in the confrontation of 
the calculated polarizabilities and the experimental re- 
sults can be considered of less importance. It is worth 
to note that the theoretical polarizabilities oscillate and 
manifest a decreasing trend such as the experimental val- 
ues do. It converts the experimental setup designed by 
Fedrigo et al. in a valuable technique to study the elec- 
tronic properties of small silver clusters. In Fig. [31 we 
have plotted the average bond length of the lowest-energy 
structures collected in Table Q] versus the cluster size. We 
have observed an abrupt change of the first-neighbor dis- 
tance from Ag2 to Ag3 and from Ag6 to Ag7, which is 
in our opinion a consequence of the structural transition 
from 1-D to 2-D and from 2-D to 3-D, respectively. The 
structural transition clearly affects the static polarizabil- 
ities as we will see below. Moreover, we can see that in 
general the average bond length approaches to the exper- 
imental value of equilibrium interatomic distance of fee 
silver solid (2.89 A) as the cluster size gets bigger [37|]. It 
is a consequence of the very important role of the surface 
effects in small clusters where most of the atoms belong 
to the surface. 

As far as the static polarizabilities displayed in Fig.[5]is 
concerned, we observe an odd-even oscillation of the cal- 
culated polarizability per atom in function of the cluster 
size ranging from the dimer up to the hexamer. This fact 
is characteristic of clusters of atoms with an odd num- 
ber of electrons and specially for atoms with a closed 
d shell and a single valence electron like in the case 
of noble- metals (Cu, Ag, and Au), or the closely re- 
lated alkali metals (Li, Na) [38|. The even-odd oscil- 
lation of the polarizability up to the hexamer is due 
to the even-odd oscillation of the HOMO-LUMO gap 
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TABLE I: Average first-neighbor distance and relative energy of Ag n cluster isomers with 2 < n < 23. The vertical ionization 
potential is compared with the data from Ref. [2f| . The geometry notation is that of Fig. [T] 



Ag 2 I > AgiQ_4 AgiQ_ 5 I > Ag23_0 

Vertical IP Vertical IP 



Cluster 


d 




Calc. 


Exp. 


|AIP| 


Cluster 




(1 


AEjjFT 


Calc. 


Exp. 


|AIP| 




(A) 


(eV) 


(eV) 


(eV) 


(eV) 




(A) 


(eV) 


(eV) 


(eV) 


(eV) 


2_0 


2.53 


0.000 


7.73 


7.60 


0.13 


10.5 


2 


.78 


0.397 


6.63 


6.25 


0.38 


3.0 


2.79 


0.000 


5.67 


6.20 


0.53 


10.6 


2 


.81 


0.495 


6.26 


6.25 


0.01 


4.0 


2.71 


0.000 


6.54 


6.65 


0.11 


10.7 


2 


.81 


0.506 


6.24 


6.25 


0.01 


4_1 


2.63 


0.242 


6.51 


6.65 


0.14 


10.8 


2 


.78 


0.553 


6.36 


6.25 


0.11 


4_2 


2.56 


0.563 


6.93 


6.65 


0.28 


10.9 


2 


.77 


0.617 


5.99 


6.25 


0.26 


5.0 


2.72 


0.000 


6.33 


6.35 


0.02 


10.10 


2 


.78 


0.626 


6.04 


6.25 


0.21 


5.1 


2.77 


0.430 


6.16 


6.35 


0.19 


10.11 


2 


.74 


0.703 


6.34 


6.25 


0.09 


6.0 


2.72 


0.000 


7.22 


7.15 


0.07 


10.12 


2 


.80 


0.720 


6.00 


6.25 


0.25 


6.1 


2.72 


0.079 


7.09 


7.15 


0.06 


10.13 


2 


.79 


1.151 


6.14 


6.25 


0.11 


6.2 


2.80 


0.534 


6.67 


7.15 


0.48 


10.14 


2 


.75 


1.133 


5.90 


6.25 


0.35 


7.0 


2.79 


0.000 


6.36 


6.40 


0.04 


11.0 


2 


.81 


0.000 


6.66 


6.30 


0.36 


7.1 


2.76 


0.170 


6.44 


6.40 


0.04 


11.1 


2 


.81 


0.268 


6.32 


6.30 


0.02 


8.0 


2.80 


0.000 


6.44 


7.10 


0.66 


12.0 


2 


.82 


0.000 


7.01 


6.50 


0.51 


8.1 


2.76 


0.006 


7.33 


7.10 


0.23 


12.1 


2 


.79 


0.117 


6.91 


6.50 


0.41 


8.2 


2.78 


0.181 


6.66 


7.10 


0.44 


12.2 


2 


.80 


0.161 


6.72 


6.50 


0.22 


9.0 


2.80 


0.000 


6.53 


6.00 


0.53 


12.3 


2 


.81 


0.411 


6.65 


6.50 


0.15 


9.1 


2.81 


0.039 


5.83 


6.00 


0.17 


12.4 


2 


.81 


0.417 


6.50 


6.50 


0.00 


9.2 


2.81 


0.039 


5.97 


6.00 


0.03 


12.5 


2 


.81 


0.422 


6.49 


6.50 


0.01 


9.3 


2.78 


0.052 


5.72 


6.00 


0.28 


12.6 


2 


.81 


0.652 


6.42 


6.50 


0.08 


9.4 


2.78 


0.056 


5.83 


6.00 


0.17 


12.7 


2 


.80 


0.802 


6.36 


6.50 


0.14 


9.5 


2.78 


0.070 


5.91 


6.00 


0.09 


12.8 


2 


.80 


1.039 


6.58 


6.50 


0.08 


9.6 


2.77 


0.089 


5.99 


6.00 


0.01 


12.9 


2 


.71 


1.039 


6.53 


6.50 


0.03 


9.7 


2.78 


0.145 


5.77 


6.00 


0.23 


13.0 


2 


.87 


0.000 


5.75 


6.34 


0.59 


9.8 


2.79 


0.154 


5.70 


6.00 


0.30 


13.1 


2 


.72 


0.305 


6.23 


6.34 


0.11 


9.9 


2.75 


0.182 


6.45 


6.00 


0.45 


13.2 


2 


.62 


3.568 


5.88 


6.34 


0.46 


9.10 


2.77 


0.204 


6.31 


6.00 


0.31 


14.0 


2 


.83 


0.000 


5.87 


6.73 


0.86 


9.11 


2.79 


0.211 


5.86 


6.00 


0.14 


15.0 


2 


.83 


0.000 


5.82 


6.40 


0.58 


9.12 


2.78 


0.214 


6.55 


6.00 


0.55 


16.0 


2 


.82 


0.000 


5.79 


6.57 


0.78 


9.13 


2.77 


0.218 


6.11 


6.00 


0.11 


17.0 


2 


.82 


0.000 


5.84 


6.45 


0.61 


9.14 


2.77 


0.221 


5.91 


6.00 


0.09 


18.0 


2 


.82 


0.000 


5.95 


6.53 


0.58 


10.0 


2.82 


0.000 


6.58 


6.25 


0.33 


19.0 


2 


.85 


0.000 


5.35 


6.20 


0.85 


10.1 


2.77 


0.080 


7.27 


6.25 


1.02 


20.0 


2 


.85 


0.000 


5.39 


6.45 


1.06 


10.2 


2.78 


0.158 


6.68 


6.25 


0.43 


21.0 


2 


.84 


0.000 


5.36 


5.90 


0.54 


10.3 


2.79 


0.205 


6.51 


6.25 


0.26 


22.0 


2 


.82 


0.000 


5.56 


6.04 


0.48 


10.4 


2.78 


0.221 


6.77 


6.25 


0.52 


23.0 


2 


.84 


0.000 


5.42 


6.03 


0.61 



(see Table [TT| in conjunction with the symmetry of the 
ground state geometries). Roughly speaking an increase 
of the HOMO-LUMO gap is on the side of a chemi- 
cal stabilization of the cluster but the chemical stabil- 
ity is favored by three-dimensional spherical structures 
or highly symmetric two-dimensional geometries which 
lead to a lost of the static polarizability. It is clearly 
reflected in Fig. 0Ja), where clusters with small HOMO- 
LUMO gaps have larger polarizability than those with 



large gaps. Thus, the static polarizability oscillates in- 
versely as HOMO-LUMO gap does. The odd-even oscil- 
lating trend is broken at n=7 because of the shape tran- 
sition from the planar to the compact three-dimensional 
structures and is reflected by a significant decrease in 
the polarization of Ag7 despite the fact that the HOMO- 
LUMO gap decrease in relation to Age and Ags. In this 
case, the symmetry of the structure dominates over the 
HOMO-LUMO gap. 



TABLE II: Calculated static response and electronic structure 
properties of the lowest-energy DFT-optimized Ag n clusters. 
The disproportionation energy is denoted by A 2 E„ and A£ 
stands for the HOMO-LUMO gap. The mean static polar- 
izability per atom a a t and the polarizability anisotropy per 
atom Aa a t were calculated under the influence of an external 
electric field of strength 0.0005 a.u.. The absolute value of 
the dipole moment is denoted by \x. 



cluster Symmetry 



A 2 E„ 
(eV) 



A£ a at 
(eV) (A 3 /atom) 



Aa a t 
(A 3 /atom) 



(D) 



1.0 






1.24 


6.85 





0.06 


2_0 


Dooft, 


0.74 


2.08 


6.88 


6.33 


0.18 


3.0 




-0.97 


0.70 


7.86 


7.21 


0.58 


4.0 


D 2h 


0.28 


0.89 


7.15 


7.83 


0.13 


5.0 


C 2 v 


-0.56 


0.55 


7.35 


6.74 


0.17 


6.0 




0.43 


2.19 


7.15 


5.69 


0.17 


7.0 


D 5h 


-0.42 


0.43 


6.47 


2.07 


0.30 


8.0 


D 2d 


0.79 


1.72 


6.32 


1.22 


0.46 


9.0 


C s 


-0.73 


0.37 


6.62 


3.49 


0.54 


10.0 


v 2d 


0.43 


0.97 


6.54 


3.99 


0.77 


11.0 


Ci 


-0.47 


0.28 


6.51 


4.05 


0.41 


12.0 


c s 


1.21 


0.83 


6.41 


3.73 


0.60 


13.0 


lh 


-0.73 


0.62 


5.96 


0.02 


0.66 


14.0 


C3v 


-0.16 


0.43 


5.75 


1.25 


0.76 


15.0 


C 2 v 


-0.02 


0.26 


5.69 


1.56 


1.20 


16.0 


c 3 


-0.05 


0.21 


5.62 


1.53 


1.18 


17.0 


C 2 


0.17 


0.20 


5.56 


1.51 


1.32 


18.0 


a 


1.11 


0.65 


5.53 


1.66 


0.94 


19.0 




-1.10 


0.12 


5.67 


2.22 


1.15 


20.0 


C 2 v 


0.33 


0.09 


5.73 


1.66 


1.16 


21.0 


Ci 


0.21 


0.15 


5.73 


1.45 


0.94 


22.0 


c s 


0.12 


0.14 


5.62 


1.80 


1.93 


23.0 


D 3h . 




0.19 


5.84 


3.06 


1.52 



From n=7 up to n=18, the static polarizability de- 
creases smoothly in accordance with the expectation that 
the polarizability per atom of a cluster is a quantity 
that decreases as the cluster becomes more compact and 
symmetric. It is borne out by the results displayed in 
Fig. 0Jb) and the inset of Fig. [4j a). With this purpose, 
we have defined the parameter 7 that characterizes the 
geometry as 



7 



5trl 

6na 2 



(4) 



where trl is the trace of the moment of inertia tensor of 
the clusters relative to the principal axis frame, n is the 
number of atoms of the cluster and a stands for the ra- 
dius of the silver atom (a« 1.45 A). In Fig.QJb) is shown 
that the larger clusters tend to have higher 7 because 
they are structurally more elongated and consequently 
it is consistent with the classical picture that the more 
spherically symmetric the cluster is, the less polarizable 
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FIG. 2: (Color online). Mean static polarizabilities per atom 
(open squares) and the polarizability anisotropies per atom 
(solid squares) of Ag„ clusters calculated with the FF method 
as a function of the cluster size. The solid circles represent the 
experimental measurements of the polarizabilities per atom 
taken from Ref. (35|. The solid line represents the prediction 
from the jellium model. The fitted parameters are given in 
the text. 




t — 1 — 1 — 1 — 1 — 1 — r 
15 17 19 21 23 



7 9 11 13 

Cluster size (n) 

FIG. 3: (Color online). First- neighbor distance plotted 
against the cluster size. The solid line represents the ex- 
perimental silver dimer bond length and the dashed line 
symbolizes the experimental value of the equilibrium inter- 
atomic distance of fee silver bulk. The three different areas 
stand for the one-dimensional (1-D), bidimensional (2-D) and 
three-dimensional (3-D) lowest-energy structures of small sil- 
ver cluster, respectively. 



is. An exception is found for Ag7 and Agi3 but the rea- 
son will be commented below in the case of Agi3 because 
for Ag7 the structure determines the reduction of the po- 
larizability as stated above. Despite the fact that the 
geometry is a fundamental parameter to describe the po- 
larizability, however it is not the only one. It is necessary 
to take into account the influence of the HOMO-LUMO 
gap. Generally, the polarizability is a result of the com- 
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FIG. 4: (a) The calculated HOMO-LUMO gap A£ and (b) 
the normalized moment of inertia 7 plotted against the mean 
static polarizability. The inset in the upper panel shows the 
evolution of the normalized moment of inertia versus the clus- 
ter size. In the bottom panel, the inset compares the polariz- 
ability anisotropy with 7. 



petition between the former and the latter contributions, 
as is shown in Fig. 2J Two structures (Ag 8 and Agi 3 ) 
manifest a significant reduction of the polarizability that 
is clearly reflected in the polarizability anisotropy which 
measures the symmetry or more specifically the defor- 
mation of the charge distribution under the influence of 
an external electric field in such a way that the less the 
polarizability anisotropy is, the more spherically symmet- 
ric the charge distribution is (see the inset of Fig.[4jb)). 
Thus, the polarizability anisotropy for Ags and Agi3 is 
clearly reduced since that Agg is a closed-shell cluster 
with a large HOMO-LUMO gap and Agi3 condensates in 



a highly symmetric structure, i.e. the icosahedral struc- 
ture (Jh). The transition from Ag§ to Agg is accompanied 
by an enhancement of the polarizability. It is caused by 
the level structure since that whenever a new level starts 
to fill, the large spatial extent of the new wave function 
contributes to the enhancement of the polarizability [3!| . 
That is the reason for the reduced value of the polar- 
izability in the case of the closed-shell structures with 
n=2,8, and 18. 

From n=19 up to n=23, the mean static polarizability 
exhibits a trend with a small positive slope. It has been 
already found in Ref . [35| . Fedrigo et al. speculates that 
this tendency is due to the ever increasing role of the 
d electrons as the cluster size grows. They state that a 
shift to red of the plasmon resonance due to d interband 
transitions corresponds to an enhancement of the polariz- 
ability. Our DFT calculations confirm this argument and 
the very important role of the d electrons as the cluster 
increases in size. In Fig. [5j we show the evolution of the 
partial density of states with the cluster size. Both d and 
sp levels gradually broaden and overlap with each other 
approximating to an electronic band in the bulk limit. 
With the aim of clarify the important role of d electrons 
with the size evolution, we have defined the energy sep- 
aration A sc i as the difference between highest occupied 
molecular orbitals belonging to 4d states and the lowest 
occupied molecular orbitals from 5s states. The energy 
separation decreases rapidly from 1.94 eV for Ag2 to 0.31 
eV for Ag6. The structural transition from planar ge- 
ometry to a three-dimensional structure gives rise to an 
increase of the value of A S( j up to 2.15 eV for Ag7. After 
that, it decreases very rapidly up to 0.07 eV for Agi 4 . For 
clusters ranging in size from n=15 up to n=23, the influ- 
ence of 4d level is so important that merges into the 5s 
state. Moreover, the small HOMO-LUMO gap collected 
in Table [TTI for n=19-23 compared to smaller clusters fa- 
vors the slightly increase of the mean static polarizability 
per atom as was commented above. 

It is expected that the jellium model can be applied 
to silver clusters because the electronic configuration of 
silver is similar to the alkali clusters, where the jellium 
model has been successfully applied. The solid line in 
Fig. [2] symbolizes the polarizability predicted by the jel- 
lium model. We have fitted our calculated mean static 
polarizabilities per atom a at to an expression given by 

m 



Uat 



in 



1/3, 



8f 



(5) 



where r ws is the Wigner-Seitz radius and 6 represents the 
spillout of the electrons from the surface of a metallic 
sphere. The values of the parameters resulting from the 
fitting to Eq. (© are r ws w 1.63 A and S ps 0.45 A which 
are close to the values r ws sw 1.58 A and S ~ 0.79 A re- 
ported in Ref. [4l[ and Ref. [l|, respectively. The bulk 
limit of Eq. ([5]) predicts a value for the bulk atomic polar- 
izability of 4.33 A 3 /atom that is lesser than those of the 
the alkali metals like for example Na which is around 9 
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Energy (eV) 



FIG. 5: (Color online). Evolution of the partial density of 
states of the the higher-lying occupied and lower-lying unoc- 
cupied levels with the cluster size. The solid, dashed, and 
dotted lines represent the contribution of the 5s, 5p, and 4d 
orbitals to the total density of states, respectively. The dotted 
vertical line represents the Fermi level. 



larizabilities from the predictions of Eq. ([5]) are "true" 
quantum effects. As commented above, it is in part due 
to the shell effects. 

In Table HH we have collected the numerical values of 
the disproportionation energy, that is defined as 



A2-E n — E n+ i + E n -x — 2E n 



(6) 



where E n is the total energy provided by our DFT cal- 
culations of the cluster with n atoms. It represents the 
relative stability of a cluster with n atoms in comparison 
to clusters with n+1 and n-1 atoms and consequently 
a peak in A2E n indicates that the cluster with size n is 
very stable because a shell has been filled. The dispropor- 
tionation energy shows that clusters with n=2,8, and 18 
have particularly stable configurations and consequently 
the polarizability is considerably reduced, as is shown in 
Fig. [2] Whenever a shell starts to fill, the polarizability 
increases and deviates from the jellium model. As the 
size of the silver clusters increases, the HOMO-LUMO 
gap becomes smaller (see Table |TTJ) so that the shell ef- 
fects are less important and the deviation of the jellium 
model is negligible. 



V. SUMMARY 

In this article the structural stability along with the 
static response properties of silver clusters in the size 
range 1 < n < 23 have been studied by means of the finite 
field method implemented in the Kohn-Sham density- 
functional methodology 12]. The IPs reported in this 
article for the lowest-energy structures are in general in 
a relatively good agreement with the experimental mea- 
surements and most of the structures predicted in this 
article as the fundamental ones were already reported 
in former publications. Likewise, the calculated polar- 
izabilities are in good agreement with the experimental 
measurements reported in Ref. 35]. The competition 
between the HOMO-LUMO gap and the structural sym- 
metry on one side or the shell structure and the dispro- 
portionation energy on the other side are the quantum- 
mechanical effects that deviates the calculated polariz- 
abilitics from the jellium model. For bigger cluster sizes 
the quantum-mechanical effects can be considered of less 
importance, and therefore both theoretical approaches, 
i.e. the ab initio DFT calculations and the jellium model 
approach each other. 



A 3 /atom [42j]. As commented above, it is due to the ever 
increasing role of the d electrons since that the screening 
of d electrons (core polarization) tends to reduce the po- 
larizability. Despite that the jellium model in the spillout 
approximation predicts in average the trend of the polar- 
izability per atom in function of the cluster size, however 
it can not account for the more interesting quantum me- 
chanical effects. Thus, deviations of the calculated po- 



APPENDIX: HELLMANN-FEYNMAN 
THEOREM 

The purpose of this appendix is to show that the finite 
basis set GGA calculation holds the Hellmann-Feynman 
theorem when fully converged in the framework of den- 
sity functional theory [43| • We have restricted ourselves 
for brevity of the formulas to wave functions without spin 
polarization, however this is not a substantial restriction 
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and the extension to spin-unrestricted orbitals is straight- 
forward. 

We have selected the ansatz in which the Kohn-Sham 
orbitals are represented by linear combinations of 

atomic Gaussian-type orbitals Xji 1 ')- Thus, the orthonor- 
mal Kohn-Sham orbitals are given by: 

^i(r) =^c ii x J -(r) (A.l) 

j 

where are the corresponding molecular orbital coeffi- 
cients. With this expansion we find the following relation 
for the electronic density: 

hj 

where Pij represents an element of the density matrix, de- 
fined as = 2 J2T C °ikCjk- Using the Eq. (|A.1|) for the 
LCGTO expansions of the Kohn-Sham orbitals subject to 
the orthonormalization condition and the electronic den- 
sity described in Eq. (|A.2|) . the variationally minimized 
Kohn-Sham SCF energy expression may be written after 
some manipulation as: 

escf (a) = E P ^ H v ( A ) + \ E p v p u (y l fc 0(A) 

ij ijkl 

- 2 efe (cjfej Cjk Sjj(X) - 1) 

ijk 

+E xc (p(\),Vp(\)) (A.3) 

where Hij represents the matrix elements of the 
core Hamiltonian and they are built from the ki- 
netic and electron-nuclear interaction energies. The 
second term represents the Coulomb repulsion energy 
of the electrons and the term E xc is the XC en- 
ergy in the GGA. We use the notation (ij\kl) = 
J f ipi(l)ip ] (l)(l/ri2)ipk(2)ipi(2)dr 1 dr2, and A being any 
parameter at all which affects the Hamiltonian of the sys- 
tem. The quantities are one-electron eigenvalues for 
the occupied orbitals and S is the overlap matrix defined 

as Sij = (Xi\Xj)- 

Assuming a gradient-corrected form for the XC energy 
E xc — j .9(p, | Vp)| 2 )<ir and upon differentiation of the 



energy with respect to A we find 
V X E = E^jVA-ffij ■ + \^2Pi j Pki^x(ij\kl) 

ij ijkl 

-2^e k CHCj k V \S lj + y^^P ij ((V\i\E xc \j) 

ijk ij 

+{i\E xc \V x j) + (i\j-V\p\j) 
+2(i\-^Vx(\Vp\)\j)) (A.4) 

The XC contribution to the derivative of the energy in- 
volves derivatives of the wave function either explicitly 
or implicitly throughout the electronic density. Denoting 
a, as a parameter of the wave functions that can be an 
exponent or positions of the bases functions, the func- 
tional derivative of the wave function can be written as 
Va«) = \di/dai)(dcii/\). Thus, we can optimize all pa- 
rameters a, so that the derivative of the XC energy can 
be neglected as well as the two-electron contribution of 
the Hamiltonian operator of Eq. (| A. 3[) . Considering that 
the overlap matrix is independent of the perturbation A 
like for example in the case of an uniform external elec- 
tric field, the third term depending on the derivative of 
the overlap matrix vanishes. Consequently, 

\7 X E = (V \H) (A.5) 

which means that the fully self-consistent finite basis set 
solutions satisfy the Hellmann-Feynman theorem in the 
framework of DFT when the XC energy is approximated 
by the GGA implementation. 
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